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Heat transport has large potentialities to unveil new 
physics in mesoscopic systems. A striking illustration is 
the integer quantum Hall regime 1 , where the robustness 
of Hall currents limits information accessible from charge 
transport—. Consequently, the gapless edge excitations are 
incompletely understood. The effective edge states the- 
ory describes them as prototypal one-dimensional chiral 
fermions 3,4 - a simple picture that explains a large body 
of observations 5 and calls for quantum information exper- 
iments with quantum point contacts in the role of beam 
splitter s 6 ! 7 ' 8 ' 9 ' 10 ! 11 . However, it is in ostensible disagreement 
with the prevailing theoretical framework that predicts, in 
most situations 12 , additional gapless edge modes 13 . Here, 
we present a setup which gives access to the energy distri- 
bution, and consequently to the energy current, in an edge 
channel brought out-of-equilibrium. This provides a strin- 
gent test of whether the additional states capture part of 
the injected energy. Our results show it is not the case and 
thereby demonstrate regarding energy transport, the quan- 
tum optics analogy of quantum point contacts and beam 
splitters. Beyond the quantum Hall regime, this novel spec- 
troscopy technique opens a new window for heat transport 
and out-of-equilibrium experiments. 

The integer quantum Hall effect, discovered nearly 
thirty years agoi, has recently experienced a strong 
revival driven by milestone experiments toward quan- 
tum information with edge states&2^. Beyond 
Hall currents, new phenomena have emerged that 
were unexpected within the free one-dimensional chi- 
ral fermions (1DCF) model. The on-going debate 
triggered by electronic Mach-Zehnder interferometers 
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vividly illustrates the gaps in our 



understanding. Coulomb interaction is seen as the key 
ingredient. In addition to its most striking repercus- 
sion, the fractional quantum Hall effect, the edge re- 
construction (ER) turns out to have deep implications 
on edge excitations. This phenomenon results from the 
competition between Coulomb interaction that tends to 
spread the electronic fluid, and the confinement poten- 
tial: as the latter gets smoother, the non-interacting edge 
becomes unstable* 19 . Theory predicts new branches of 
gapless electronic excitations in reconstructed edge a 13 ' 20 , 
which breaks the mapping of an edge channel (EC) onto 
1DCF and, possibly, the promising quantum optics anal- 
ogy. For most edges realized in semi-conductor hetero- 
junctions (except by cleaved edge overgrowth 2 ^) , ER re- 
sults in wide compressible ECs separated by narrow in- 
compressible strips 12 and the new excited states are over- 
all neutral internal charge oscillations across the ECs 
widtbM. 

In practice, the predicted additional neutral modes are 
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Figure 1: Experimental implementation of non- 
equilibrium edge channel spectroscopy, a, Schematic 
description of the energy distributions fo,s{E) spectroscopy 
with a single active electronic level of tunable energy Ei ev (VG) 
in the quantum dot (QD). b, The current Tqd {dlQu/dVc) is 
proportional to f s {E)-f D {E) (d(f s {E)-f D (E))/dE) ignor- 
ing variations in tunnel rates and tunneling density of states, 
c, E-beam micrograph of the sample. Surface metal gates 
appear brighter. Electronic excitations propagate counter 
clockwise along two edge channels (EC) of the quantum Hall 
regime. The outer EC (continuous white lines) is partly trans- 
mitted (dashed lines) across quantum point contact (QPC) 
and QD. The inner EC (not shown) is always reflected. The 
QPC is used to drive out-of-equilibrium the drain outer EC. 
Gates partly covered by the insets are grounded and do not 
influence the electron paths. Left inset: non-interacting elec- 
trons prediction for /d(-E) in the outer EC at output of QPC. 



transparent to Hall currents. More surprisingly, a lin- 
ear I-V characteristic is frequently observed for tunnel 
contacts (different behaviors were also reported, e.gj^ 2 .) 
whereas a non-linear characteristic is predicte d 23 ! 24 ' 25 . 
This contradiction is resolved by assuming ad-hoc that 
only rigid displacements of compressible ECs are excited 
by tunnel events, and not internal excitation a 13 ' 23 ! 24 ' 26 . 
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The rigid displacement model arguably relies on the over- 
riding strength of Coulomb interaction that tends to or- 
thogonalize bare tunneling electrons and correlated elec- 
tronic fluids^. However, the above argument does not 
hold at arbitrary transmission probabilities, where mul- 
tiple electrons processes occur. Therefore, the role of 
predicted internal excitations has to be determined ex- 
perimentally. The present work provides such a test. An 
EC is driven out-of-equilibrium with a quantum point 
contact (QPC) of arbitrary transmission, possibly excit- 
ing internal modes. A short distance away, the result- 
ing energy distribution f(E) is measured with a tunnel 
coupled quantum dot (QD) expected to probe only rigid 
displacement excitations, hereafter called quasiparticlcs. 
Consequently, the amount of energy injected into inter- 
nal modes at the QPC would appear as an energy loss in 
f(E). 

Measurements of the energy distribution in mesoscopic 
devices were first performed in 1997 on metallic circuits 
using a superconducting tunnel probe^l. More recently, 
the regime of high magnetic fields was accessed using dy- 
namical Coulomb blockade 2 ^. In two-dimensional elec- 
tron gas systems, non-Fermi energy distributions could 
not be measured since transferring the techniques de- 
veloped for metal circuits is technically challenging (al- 
though hot electrons have been detected, e.gi 2 ^). Regard- 
ing the quantum Hall regime (QHR), state of the art is 
the very recent qualitative probe of heating 3 ^. Here, we 
demonstrate that f(E) can be fully extracted from the 
tunnel current across a QD. In the sequential tunneling 
regime, the discrete electronic levels in a QD behave as 
energy filter o 31 ' 32 , as previously demonstrated with dou- 
ble QDs 33 . Assuming a single active QD level of energy 
Ei ev , and ignoring the energy dependence of tunnel rates 
and tunneling density of states in the electrodes, the QD 
current reads 



^QD = lQB X (fs(Ei e 



fD(Elev)), 



(1) 



where the subscript S (D) refers to the source (drain) 
electrode; fs,D are the corresponding energy distribu- 
tions and Iqp x is the maximum QD current. In prac- 
tice, fs.D are obtained separately by applying a large 
enough source-drain voltage (Fig. [T^,b) and the probed 
energy Ei ev — Eq — erjcVc is swept using a capac- 
itively coupled gate biased at Vq, with rjo the gate 
voltage-to-energy lever arm and Eq an offset. Raw data 
9Iqd/9Vg measured by lock-in techniques are propor- 
tionalto df D ,s(E)/dE. 

A tunable non-Fermi energy distribution is generated 
in an EC with a voltage biased QPC. Similar setups 
were used previously to create imbalanced electron pop- 
ulations between co-propagating ECs^i, each character- 
ized by a cold Fermi distribution. Only in a very re- 
cent experiment 3 -^, was an EC heated up. Beyond heat- 
ing, f{E) is here controllably tuned out of equilibrium. 
Let us consider one EC and assume it can be mapped 
onto non-interacting 1DCF. According to the scattering 
approach 3 -^, the energy distribution at the output of a 
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Figure 2: Equilibrium edge channel spectroscopy and 
quantum dot characterization. Temperatures Tf it are 
obtained by fitting the QD-drain dip (o) and QD-source peak 
(■) in <97qd/<9Vg (resp. left and right peaks in bottom right 
inset), and from the single OIqd /8Vd peak at Vd = (*) (raw 
data not shown) with r/c = 0.057. We assumed a single active 
level in the QD and Fermi energy distributions in source and 
drain. Errors within T± 10% are enclosed between the dashed 
lines. Top left inset: Surface plot of <9Jqd/<?Vg (negative is 
brighter, positive darker) measured at T = 30 mK for an 
outer (inner) drain EC biased at Vd (—88 /tV). The intense 
stripes slopes yield r\G = 0.052 ± 9%. Enclosed dashed lines 
outline the small contributions of other electronic levels. 



QPC of transmission t is a tunable double step (Fig. [IF,, 
left inset) 



fD(E)=Tf m (E) + (l-T)f D2 (E), 



(2) 



where foi (/ds) is the equilibrium Fermi distribution 
function in the partially transmitted (reflected) incoming 
EC of electrochemical potential shifted by eVoi (eVo2)- 
In presence of ER, the above energy distribution applies 
to the quasiparticles if internal modes are not excited at 
the QPC. Else, there are no theoretical predictions be- 
cause a QPC is very difficult to treat non-perturbatively 
in their natural bosonic formalism. 

The sample shown in Fig. [TJ: was tailored in a two- 
dimensional electron gas realized in a GaAs/Ga(Al)As 
heterojunction, set to filling factor 2 and measured in 
a dilution refrigerator of base temperature 30 mK. The 
experiment detailed here focuses on the outer EC rep- 
resented as a white line. The inner EC (not shown) 
is fully reflected by the QPC and the QD. We checked 
that charge tunneling between ECs is negligible along the 
0.8 /mi propagation length from QPC to QD. 

We first perform a standard non-linear QD 
characterization 32 (Fig. [2J top left inset). The two 
large signal stripes are frontiers of consecutive Coulomb 
diamonds and are accounted for by a single active QD 
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Figure 3: Spectroscopy of an edge channel tuned out- 
of-equilibrium with the conductance of a QPC. The 

QPC's voltage bias is here fixed to SVd = Vdi — Vd% = 36 yN . 
Data shown as (+) in (c) and (d) were obtained in a second 
cooldown with a renewed rjc = 0.059 ± 7%. a, Measured 
Gqpc versus Vqpc applied to the lower left metal gate in 
Fig. [Tp. Symbols are data points obtained during the corre- 
sponding EC spectroscopy. The continuous line was measured 
earlier with SVd = 0. b, Measured dlqu/dVc (symbols) for 
the QD-drain contribution (data have been aligned in Vg and 
normalized to Iq£ x )- The two dips correspond to a double 
step energy distribution fr>(E). c, Symbols are Tfu obtained 
from fits of the data (continuous lines in (b)) assuming fo(E) 
is the weighted sum of two Fermi functions. We find an ac- 
curate agreement with the non-interacting ID chiral fermions 
(1DCF) model prediction Tfu = Gqpch/e 2 . d, Generalized 
non-equilibrium temperature (symbols, see text) extracted 
from the data and theoretical prediction for free 1DCF (con- 
tinuous line). Error bars are dominated by uncertainties on 
r\G- e, Energy distributions obtained by integrating the data 
in (b). 



level. Small contributions of three additional levels 
of relative energies {-95,30,130} fieV are also visible. 
The lever arm extracted from the stripes' slopes is 
r, G ~ 0.052 ± 9%. 

Then, we test the spectroscopy with known Fermi func- 
tions by measuring <9/qd/<9Vg(Vg) at Vdi — Vdi — 
—88 fiV for several temperatures (Fig. [21 bottom right 
inset). By fitting these data with Eq.[T]using Fermi func- 
tions, we extract a fit temperature scaled by the lever 
arm Tfa/r/G- The value t\g — 0.057, compatible with 
the non-linear QD characterization, is found to repro- 
duce best the mixing chamber temperature T with Tfu- 
The drain (o) and source (■) fit temperatures are shown 
in Fig. [^together with Tfu obtained using the standard 
procedure^ from 9/q D /9Vd(Vg) at Vd — (*). We find 



deviations mostly within ±10% (dashed lines in Fig. [2]) 
except for a saturation at Tfu ~ 50 mK possibly due 
to a higher electronic temperature. In the following, we 
use r\Q — 0.057 obtained here in the same experimental 
configuration as to measure unknown f(E)s. 

Electrons are now driven out-of-equilibrium in the 
drain outer EC. In the following, the electrode D2 and 
the inner drain EC are voltage biased at Vdi — —88 /iV 
and the source ECs are emitted by a cold ground. 

First, the bias voltage across the QPC is set to 5Vd = 
Vd\ — Vd2 = 36 fiV and its conductance Gqpc = re 2 /h is 
tuned by applying Vqpc to the bottom left gate in Fig.[Tp 
(see Fig. [3£i). Note that at 30 mK, we find the transmis- 
sion t is constant within 2% with the QPC voltage bias 
below 36 /iV. Typical sweeps 9/qd /9Vg(Vg) and the cor- 
responding Jd{E) are shown in Figs.[3p and [3p, respec- 
tively. The QD-drain negative contribution transforms 
from a single dip at t = {0,1} into two dips separated 
by a fixed gate voltage and whose relative weights evolve 
monotonously with r e]0,l[. Continuous lines are fits 
with Eq. fusing for Jd\.D2 two Fermi functions shifted 
by a fixed energy and weighted by the factors Tfu and 
1 — Tf a . The values of Tfu are found to deviate by less 
than 0.03 from the measured transmission t (Fig. [3p), in 
accurate agreement with the free 1DCF model. Symbols 
(+) in Fig. [3] correspond to data obtained in a second 
cooldown. 

In a second step, the QPC transmission is fixed to 
t w 0.5 and the bias voltage SVd is changed. Typical 
raw data are shown in Fig. 2^. These were obtained 
in a third cooldown with a QD renewed by the ther- 
mal cycle showing no signs of additional QD levels in 
the probed energy range. The single dip in the QD- 
drain contribution (bright) at SVd = splits in two 
similar dips that are separated by a gate voltage differ- 
ence proportional to SVd- Meanwhile, the QD-source 
peak (dark) is mostly unchanged but slowly drifts par- 
allel to one QD-drain dip due to the capacitive coupling 
between drain and QD. In the first cooldown, Vdi was 
kept within [—106,-34] /iV to minimize complications 
related to additional QD levels (lower bound) and to en- 
sure well-separated source and drain contributions (up- 
per bound). Symbols in Figs. Hp and [3£ are, respec- 
tively, data and extracted Jd{E) for the QD-drain con- 
tribution at 5V D = {-18,0,18,27,36,45,54} (iV and 
t = 0.58. Continuous lines in Fig.QJ) are fits with Eq. [5] 
using the measured r and for /di,D2 two Fermi functions 
shifted in energy by the fit parameter —eriGSVc- The 
resulting tjgSVg are plotted as symbols versus SVd in 
Fig. [3p. Those obtained in the third cooldown are shown 
as (*) using the renewed lever arm t\g = 0.062. We 
find tjqSVg — SVd as expected from the non- interacting 
1DCF model. Deviations are always smaller than 8 fiV 
(5 fiV) for the first (third) cooldown, a reasonable agree- 
ment regarding uncertainties in t]g of ±10% (±5%). 

In the two experiments above, we found the measured 
quasiparticle f(E)s verify predictions of the scattering 
approach. In order to establish the analogy QPC-beam 
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Figure 4: Spectroscopy of an edge channel tuned out- 
of-equilibrium with the voltage across a QPC. The 

QPC's conductance is here fixed to Gqpc ~ 0.5e 2 /h. Data 
of (a) and (*) in (c) and (d) were obtained in a third 
cooldown with a renewed no = 0.062 ±5%. a, Surface plot of 
dlQD / 8Vg{5Vd, Vg)- The QD-drain contribution (brighter) 
splits into two similar dips separated by a gate voltage dif- 
ference proportional to SVd- b, Measured dlQv/dVG (sym- 
bols) for the QD-drain contribution. Data have been shifted 
vertically for clarity, and horizontally to align the peak corre- 
sponding to the fixed potential Vd2- Continuous lines are fits 
assuming Jd(E) is the weighted sum of two Fermi functions, 
c, Symbols are the fit parameters ncSVa- The continuous 
line is the prediction for non-interacting 1DCF. d, General- 
ized non-equilibrium temperature (symbols) and theoretical 
prediction for free 1DCF (continuous line). The good agree- 
ment data-prediction demonstrates that internal modes are 
not excited at the QPC within our experimental accuracy, e, 
Energy distributions obtained by integrating the data in (b). 

splitter one also needs to demonstrate that internal EC 
modes are not excited. A direct test consists in extract- 
ing the quasiparticlc heat current J|? from the data, and 
comparing it with the full edge excitations heat current 
Je obtained from power balance considerations (see Sup- 
plementary Information for details): 

MT = 0) = ^r(l-r). (3) 



The cancellation vv — l/h of velocity (v) and density of 
states per unit length and energy (v) that applies to the 
1DCF quasiparticles permits us to obtain Jjf from the 
measured f(E) without any sample specific parameters: 



J E=lJ( E - - 0(1* - E))dE, (4) 



with [i the electrochemical potential and 0(E) the step 
function. Consequently, we measure quantitatively the 
quasiparticle heat current. The result of this pro- 
cedure is shown as symbols in Figs. and [JJl us- 
ing the generalized non-equilibrium temperature T qp = 
y QhJW /^kg together with the prediction Tidcf — 
y/T 2 + t(1 — T)3(e5Vo l^ks) 2 if none of the injected 
power is carried on by internal modes (continuous lines). 
We find a good agreement T qp ~ Tidcf without fit- 
ting parameters and essentially in or close to error bars. 
Hence, within our experimental accuracy, the propaga- 
tive internal modes do not contribute to heat transport 
and therefore are not excited. Note that the relatively 
small observed deviations are cooldown dependent, which 
points out the QD. Indeed, the data can be more accu- 
rately accounted for including a second active QD level 
(see Supplementary Information). Last, preliminary data 
show a significant energy redistribution with the inner 
EC for propagations longer than 2 //m in the probed 
energy range. Therefore, the observed small discrepan- 
cies could also result from the finite 0.8 /zm propagation 
length. 

Overall, we demonstrate that QPCs in the QHR 
are tunable electrical beam splitters for one-dimensional 
fermions, i.e. rigid edge channel displacements, i) by 
comparing the energy distribution at a QPC output with 
predictions of the scattering approach 5 , and ii) by show- 
ing that internal edge channel modes are not excited. 
This does not only rule out non-ideal QPC behaviors 
to explain the surprising phenomena observed on elec- 
tronic Mach-Zchnde r 8 i 15 ' 16 i 17 . It also establishes a solid 
ground for future quantum information applications with 
edge states. Finally, an essential part of this work is the 
demonstration of a new technique to measure the fun- 
damental energy distribution function. It makes f(E) 
accessible for most systems were quantum dots can be re- 
alized. We expect it will trigger many new experiments 
dealing with heat transport, out-of-equilibrium physics 
and quantum decoherence. 
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I. EXPANDED DISCUSSION ON ENERGY AND 
HEAT TRANSPORT BY EDGE EXCITATIONS 



In this section, after a schematic summary of the exper- 
imental principle, we first detail how the energy within 
the probed electronic excitations can be extracted from 
the energy distribution function f(E) and we define a 
generalized temperature T qp for non-equilibrium situa- 
tions. Second, we detail the link between energy within a 
system of ID chiral fermions, such as the 'quasiparticlcs' 
(i.e. the rigid edge channel displacement excitations), 
and energy transport. In this case, distribution function 
and heat current are directly related, without any sam- 
ple specific parameters. Third, we extract from general 
power balance considerations the total excess energy cur- 
rent carried by all edge excitations. It is compared to the 
energy current carried only by the probed quasiparticles. 
This permits us to rule out any contribution of internal 
modes to the energy transport within our experimental 
accuracy. 



A. Schematic summary of experiment principle 

Figure[T]recapitulates the logical architecture that per- 
mits us to rule out the excitation of internal EC states 
and to show that QPCs obey the scattering model for 
quasiparticles 1 . This demonstrates the quantum optics 
analogy, regarding energy transport, between a QPC and 
a tunable beam splitter for quasiparticles. 

The hypothesis HI relies on the linear I-V character- 
istic for tunnel contacts that is mostly observed (includ- 
ing here) in the integer QHR£i2iii£i£. The hypothesis H2 
holds for each edge channel's quasiparticle branch ac- 
cording to^^i 7 -. Note that the hydrodynamic approach 
considers specifically rigid EC displacements 7 . The hy- 
pothesis H3 is a prediction for internal EC modes^^. 



B. Energy density of non-equilibrium fermions 

At equilibrium and temperature T, the energy distri- 
bution function f(E) is a Fermi function and the typ- 
ical energy of electronic excitations is ksT. In non- 
equilibrium situations, the distribution function is in gen- 
eral not a Fermi function. Nevertheless, the energy den- 
sity E qp within a system of fermionic quasiparticles can 



Working hypothesis 

H1: Tunnel contacts do not couple to internal EC modes 

H2: "Quasiparticles" can be described as 1D chiral fermions 
(quasiparticles = rigid EC displacement excitations) 

H3: Internal EC modes are propagative 

Implications of hypothesis 

H1 => 11: f(E) measured with a tunnel coupled QD is that 
of quasiparticles 

H2 + 11 => 12: Quasiparticle heat current J E P given by Eq. 4 

H3 => 13: If 3 excited internal states at QPC then this implies 
excess J E P < excess edge excitations heat current J E 

Experimental observations 

E1: Measured f(E) complies with Eq. 1 
E2 (using 12): Excess J E P ~ excess J E 
Conclusions 

E1 + 11 => C1: Quasiparticles f(E) obey the scattering model 
of QPCs 

E2 + 13 => C2: Internal modes are not excited at the QPC 

C1 + C2 => C3: Regarding energy transport, QPCs in int. QHR 
are electrical analogues for quasiparticles of 
tunable beam splitters 

Figure 1: Schematic summary of experiment principle. 



be obtained from their energy distribution: 
E qp =u f (n-E)(l-f(E))dE 

I' oc 

+ u (E- ti)f{E)dE, 
J/j, 

=v[(E- »)5f(E)dE, (1) 

with v the density of states per unit length and energy, 
here assumed constant for energies near the local electro- 
chemical potential \x\ and Sf(E) = f(E) —9(fi — E), with 
0(E) the step function, corresponding to variations in the 
energy distribution relative to the filled Fermi sea. The 
second line and the right-hand side of the first line corre- 
spond to the energy density contribution of, respectively, 
electron and hole like quasiparticles. In the important 
case of a Fermi function at temperature T, the energy 
density obtained from Supplementary equation [1] is the 



standard textbook value 9 : 



E„ 
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y{k B Tf 



(2) 



In practice, it is more convenient to use the electronic 
temperature T qp generalized to non-equilibrium situa- 
tions, which is independent of the density of states, rather 
than the energy density E qp : 



~2 ' kB - 



(3) 



The energy density and generalized electronic tempera- 
ture can be computed analytically for the 'double step' 
energy distributions detailed by Eq. (2). One finds: 



L QP/ 1 ' = g 1" A 1 - J 2 : 



^ P =\ T 2 + 3t(1 - r) 



(4) 
(5) 



C. Energy transport by ID chiral fermions 

We expect from the rigid displacement model^^^ 
supported by tunneling density of states experiments, 
that the energy distribution, measured here with a 
weakly coupled quantum dot, is that of rigid displace- 
ment excitations (the quasiparticles) and not of internal 
edge channel excitations. Moreover, it was shown that 
the quasiparticles of each edge channel can be mapped 
onto a branch of ID chiral fermions^^i Therefore, the 
measured energy current J|? along an edge channel cor- 
responds to that of quasiparticles and reads: 



vE qp = ^-(k B T qp ) 2 , 



(0) 



with v the drift velocity and h the Plank constant. Note 
first that the right-hand side expression of heat current 
holds even if the density of states v depends on energy, 
as long as the very robust ID velocity-density of states 
cancelation vv = 1/h is obeyed^. Second, as expected 
the above expression obeys the Wiedemann-Franz law. 
Third, we point out that the energy flow is directly given 
by the generalized temperature T qpi without any sample 
specific parameters. Therefore, by measuring the energy 
distribution function in our experiment, we probe quan- 
titatively the quasiparticle energy current. 

Assuming the quasiparticle distribution function is 
given by Equation (2), the corresponding energy flow is: 



(>/? 



(k B T) 2 +t(1-t) 



(eSVnf 
2h 



(7) 



D. 



Excess energy transport by all edge excitations 
and comparison with that by quasiparticles 



We consider the simplified circuit at filling factor 1 (one 
edge channel) shown in Supplementary figure [21 




Figure 2: Simplified schematic circuit used for power 
balance considerations. The Landau level filling factor is 
here set to one. The edge channel is shown as a white line and 
the propagation direction of electronic excitations is indicated 
by an arrow. The dissipated power V 2 Gqpc = V 2 re 2 /h is 
injected into the source (S) and drain (D) electrodes by the 
corresponding incoming edge channels. 



Here, the total power P = V 2 Te 2 /h, with t the 
quantum point contact transmission, is dissipated in the 
source (S) and drain (D) electrodes (nearby the hot spots 
shown as red areas in Supplementary figure^ see e.g^i). 
It can be decomposed into two contributions: 



P = V 2 re 2 /h = P 5p 



P P , 



(8) 



The first one (Psfi) corresponds to the energy injected 
into the drain and source electrodes due to the electro- 
chemical potential difference d/i between electrode and 
corresponding incoming edge. The edge electrochemical 
potential is defined as that of a floating electrode inserted 
in its path, in the spirit of the 'measurement reservoir' 
model (see e.gJ^). At unity transmission r = 1, this 
'electrochemical power' is the only contribution to the 
dissipated power P$ p {t = 1) = P = (eV) 2 /h. In general, 
the electrochemical power injected by each edge in its 
output electrode is (5[i) 2 /2h. At arbitrary transmission 
t, the electrochemical potential difference at the input of 
both the source and drain electrodes is \S/i\ = re\V\ and 
one finds: 



Psn = {reVf/h. 



(9) 



The second contribution P e d ge — 2(Jg* — J™*) corre- 
sponds to the difference between the incoming J|f and 
outgoing Jjj energy current carried on by all edge ex- 
citations, respectively in and out the corresponding elec- 
trode. The factor two here accounts for the two elec- 
trodes. Note that P e dge corresponds to the amount of 
energy that would be absorbed by two floating reservoirs 
inserted along the path of the edges incoming to source 
and drain electrodes and thermalized at the same temper- 
ature as their corresponding electrode. This contribution 
vanishes at zero transmission and also at unity transmis- 
sion, as long as the drain and source electrodes are at 
the same temperature since in that case J^ 1 = Jg u *. At 
intermediate transmissions, P e dge is obtained from Sup- 
plementary equations [H] and 



edge 



P-Ps^ = r(l - r) 



(eSVn 



(10) 
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This last quantity is identical to the one obtained from 
the double step distribution function (using J l £ given 

by Supplementary equation [7] and J™* = ^(fcel 1 ) 2 ). 
Therefore the observed agreement between measured 
quasiparticle f(E) and the prediction of Eq. 2 already im- 
plies that the excess energy current is carried on by quasi- 
particles and not internal edge channel modes. However, 
the most straightforward evidence is to extract the energy 
current directly from the measured quasiparticle energy 
distribution function using Supplementary equations [IJ [3] 
and [5] The internal modes being propagative, the obser- 
vation that excess quasiparticle energy current and full 
excess edge current are similar (see Figures 3d and 4d) 
implies internal edge channel states are not excited by a 
voltage biased QPC of arbitrary transmission within our 
experimental accuracy. 



II. SUPPLEMENTARY DATA 

In this section we present supplementary data regard- 
ing cooldowns two and three, that were not shown in the 
article. 

The quantum dot calibration data are shown in Sup- 
plementary figure [3J 

Data obtained in cooldown two with a fixed bias volt- 
age SVd = —36 /iV and several values of the QPC con- 
ductance are shown in Supplementary figure [4] These 
data correspond to symbols (+) in Figure 3c, d. Note 
that here Vd2 and the inner edge channel potential are 
set to -129 /iV. 

Data obtained in cooldown three with a fixed QPC 
conductance Gqpc — 0.5e 2 //i and several values of the 
QPC voltage bias SVd are shown in Supplementary fig- 
ure As for cooldown one, Vd2 and the inner edge 
channel potential are here fixed to —88 /iV. 



III. SUPPLEMENTARY INFORMATION ON 
METHODS 

A. Measured sample and experimental techniques 

The sample was realized in a standard GaAs/Ga(Al) As 
two dimensional electron gas 105 nm below the surface, 
of density 2 10 15 m~ 2 , Fermi energy 80 K and mobil- 
ity 250 m 2 !/-^" 1 . Note that the GaAs/Ga(Al)As het- 
erojunction used here is the same one as that used for- 
merly to perform the Mach-Zchndcr experiments with 
edge states reported i n 12 i 13 . The silicon (dopant) con- 
centration in the heterojunction was adjusted to optimize 
the Hall resistance quantization. 

Conductance measurements were performed in a di- 
lution refrigerator of base temperature 30 mK. All mea- 
surement lines were filtered by commercial 7r-filters at the 
top of the cryostat. At low temperature, the lines were 
carefully filtered and thermalized by arranging them as 




T (mK) 



Figure 3: Quantum dot calibration in cooldowns two 
and three. This figure complements Figure 2 that focused on 
cooldown one. (a) and (b) are calibration data for cooldowns 
two and three, respectively. Errors in Tfu within ±7% for 
cooldown two (a) and within ±5% for cooldown three (b) 
are enclosed between the black dashed lines. Note that an 
additional QD level of relative energy —56 fieV for cooldown 
two is visible in the top left inset of (a). 



1 m long resistive twisted pairs (300 O/m) inserted in- 
side 260 fim inner diameter CuNi tubes tightly wrapped 
around a copper plate screwed to the mixing chamber. 
The sample was further protected from spurious high en- 
ergy photons by two shields, both at base temperature. 

The sample was current biased by a voltage source 
in series with a 10 MJ7 or 100 Mil polarization resis- 
tance at room temperature. Taking advantage of the well 
defined quantum Hall resistance (12.906 kS7) , currents 
across the sample were converted on-chip into voltages 
and measured with low noise room temperature voltage 
amplifiers. To limit artifacts by slowly moving charges 
nearby the QD, we systematically measured several suc- 
cessive gate voltage sweeps Iqd(Vg), checked that the 
data fall on top of each other, and verified that the sum 
rule J (<9/qd /dVc)dVG — is obeyed. 
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V Q (mV) E (neV) 

Figure 4: Spectroscopy of an edge channel tuned out- 
of-equilibrium with the conductance of a QPC in 
cooldown two. The QPC's voltage bias is here fixed to 
SVd = Vbi — Vl>2 = —36 j«V. a, Measured Gqpc versus Vqpc 
applied to the lower left metal gate in Figure lc. Symbols 
are data points obtained during the corresponding EC spec- 
troscopy. The continuous line was measured with SVd = 0. 
b, Measured OIqu/OVg (symbols) for the QD-drain contribu- 
tion (data have been aligned in Vg and normalized to Iq§ x ). 
Continuous lines are fits assuming fr> (E) is the weighted sum 
of two Fermi functions. The detailed set of used fit parameters 
is given in Supplementary table IHII c, Energy distributions 
obtained by integrating the data in (b) and using r\c = 0.059. 




Figure 5: Spectroscopy of an edge channel tuned out- 
of-equilibrium with the voltage across a QPC in 
cooldown three. The QPC's conductance is here fixed to 
Gqpc = 0.5e 2 /h. a, Measured <9/qd/<9Vg (symbols) for the 
QD-drain contribution. Data have been shifted vertically for 
clarity, and horizontally to align the peak corresponding to 
the fixed potential Vd2- Continuous lines are fits assuming 
/d(E) is the weighted sum of two Fermi functions (see Sup- 
plementary table [IV}. b, Energy distributions obtained by 
integrating the data in (a) and using r\G = 0.062. 



B. Quantum dot theoretical model 

The QD modeling follows the master equation ap- 
proach for incoherent sequential tunneling. The deriva- 
tion detailed ini^ stays valid for non-Fermi distribution 



functions. To obtain the simple expression of Equation 1 
we assumed only one QD level contributes to /qd, and 
neglected the energy dependence of the electrodes tun- 
neling density of states and of the tunnel rates in and 
out the QD. In practice, the validity of these hypothe- 
ses are checked with the non-linear QD characterization 
and by comparing mixing chamber temperatures with fit 
temperatures obtained within this framework (see Fig- 
ure 2 and Supplementary figure [3]). Note that, beyond 
sequential tunneling, a similar relationship between /qd 
and fs.D holds in the resonant transport regime but with 
distribution functions artificially smoothed by the energy 
width of the QD level. 



C. Experimental procedures 

1. Extraction of the energy distribution function 

The distribution functions f(E) are ob- 
tained by integrating with Vg the normalized 
data —(dlQn/dVG)/lQ£, x from a charge stabil- 
ity zone (Jqd = 0) on the negative side up to 
V G = -(E-E )/e VG . 



2. Extraction of the electrochemical potential 

The electrochemical potential [i enters in the evalua- 
tion of the energy density (Supplementary equation [T]) , 
and consequently in the heat current Jjf and in the gen- 
eralized quasiparticle temperature T qp . The parameter 

is here obtained very directly from the measured f(E) 
using 

H = E mm + / f(E)dE, (11) 

JEmin 

with E m i n (E max ) an energy under (above) which we 
assume f(E) = 1 (f{E) = 0). Note that /i is obtained up 
to the same unknown constant offset Eg as f(E), however 
this has no impact since only the difference E — fi plays 
a role. 



3. Estimation of uncertainties on the lever arm t\g 

Uncertainties in the gate voltage-to-energy lever arm 
conversion factor t\g are mostly responsible for our error 
bars on energy related quantities (e.g. Figs. 3d and 4c, d). 

This parameter is extracted from two independent cal- 
ibration procedures: First, we perform a non-linear char- 
acterization of the QD 15 and extract t\g from the slopes 
of the Coulomb diamond (?/qd/9Vg(Vg, Vd) (hereafter 
called procedure 1). Second, we measure <3/qd/<3Vg(Vg) 
at several temperatures and extract r\G from the scal- 
ing between fit temperatures (using Fermi functions in 
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Equation 1) and measured mixing chamber temperatures 
(hereafter called procedure 2). 

In procedure 1, uncertainties on t\q are obtained from 
the change in slopes corresponding to displacements 
equal to the full width at half maximum of the peaks 
at the highest drain-source voltage. This gives r\G — 
0.052 ± 9%, 0.062 ± 8% and 0.055 ± 9% for, respectively, 
cooldowns one, two and three. Note that these uncertain- 
ties correspond to plus and minus the full width at half 
maximum corresponding to the effective electronic tem- 
perature 50 mK (i.e. 0.05 x 3.5fcs/e ~ 15 fiV) divided 
by the maximum drain voltage Vd applied. 

In procedure 2, uncertainties in r/c are obtained from 
the dispersion in fit temperatures Tf it around measured 
temperatures T of the dilution refrigerator mixing cham- 
ber. We evaluate roughly the uncertainty by finding the 
range of t\q that permits us to account for most Tfit at 
T > 50 mK. The reader can get a direct idea of the used 
uncertainties in Figure 2 and Supplementary figure (3[ 
where the expected fit temperatures using the consid- 
ered extremal values of r/c are shown as dashed lines. 
This gives r, G = 0.057 ± 10%, 0.059 ± 7% and 0.062 ± 5% 
for, respectively, cooldowns one, two and three. 

In the article, we have chosen to use the values of rja 
and the associated uncertainties extracted using proce- 
dure 2, in the same experimental configuration as to mea- 
sure unknown f(E)s. 



4- Estimation of error bars 

The error bars in the figures take into account the finite 
signal to noise and reproducibility (i.e. the standard error 
in the average value), and, if it applies, uncertainties in 
the lever arm. 

Regarding the finite signal to noise and reproducibil- 
ity our approach is very straightforward. We perform a 
statistical analysis on the considered quantity extracted 
from up to a hundred different Vq sweeps in the ex- 
act same experimental configuration. The corresponding 
standard error is then plus/minus the mean deviation 
per sweep divided by the square root of the number of 
sweeps. This is the only contribution for the parameter 
Tfu whose typical error bars are found to be about 0.02, 
small compared to the symbol size in Fig. 3c. Note that 
in practice, we acquired the large number of sweeps nec- 
essary for an accurate statistical analysis in only a few 
realizations per experiment. In other realizations of a 
given experiment (i.e. when changing only the QPC con- 
ductance or the applied voltage bias), we assumed that 
the observed mean deviation per sweep is unchanged and 
estimate the standard error using the corresponding num- 
ber of sweeps (generally more than five). 

It turns out that, except for the immune Tfit, error 
bars are mostly dominated by uncertainties on the lever 
arm rja and therefore are proportional to the overall en- 
ergy (namely, T qp or SVd)- Indeed, the standard error on 
T qp was always found smaller than 1 mK for cooldowns 



Gqpc {e'/h) 


Tfit 


T D1 {K) 


T D2 (K) 








0.042 


0.042 


0.15 


0.18 


0.158 


0.042 


0.26 


0.25 


0.078 


0.046 


0.35 


0.33 


0.071 


0.049 


0.43 


0.42 


0.067 


0.050 


0.51 


0.48 


0.062 


0.054 


0.59 


0.58 


0.064 


0.054 


0.59 


0.61 


0.067 


0.056 


0.65 


0.65 


0.062 


0.056 


0.72 


0.73 


0.061 


0.053 


0.78 


0.77 


0.056 


0.059 


0.88 


0.86 


0.053 


0.073 


1 


1 


0.045 


0.045 



Table I: Summary of parameters used to fit the data 
shown in Figure 3b (cooldown one). The conductance 
Gqpc is measured. Note that near zero and full transmission, 
the fit temperature of the small corresponding peak is not 
very significant. In order to focus on the fit parameter r/ it , 
we chose to fix r/cSVc = 30 /iV. If tigSVg is set free, we find 
values within 30 ± 1 /iV, except at Gqpc = 0.9e 2 /h where 
r/cSVG = 39 /iV. Note that Tf it is not affected more than 
±0.03 by whether tigSVg is fixed or free. 



SV D (pV) 


T D1 (K) 


T D2 (K) 


no5V G (nV) 


-18 


0.049 


0.049 


-12 





0.040 


0.040 





18 


0.048 


0.045 


13 


27 


0.054 


0.054 


20 


36 


0.061 


0.056 


29 


45 


0.074 


0.063 


38 


54 


0.076 


0.073 


46 



Table II: Summary of parameters used to fit the data 
shown in Figure 4b (cooldown one). Here SVd is the 
applied voltage bias and the parameter t/u is set to the mea- 
sured G Q Pch/e 2 = 0.58. 



one and three, and 4 mK for cooldown two. Nevertheless, 
the full error bars shown in figures include both contri- 
butions, taken as independent from each other. 



5. Fit procedures 

We fitted the measured (9/qd / OVg)/ Tqff using Equa- 
tions 1 and 2, with Fermi functions for fr>\ and fo2- The 
full set of fit parameters is not only constituted of Tfu 
and t]gSVg shown in Figures 3c and 4c, respectively. It 
also includes the temperatures Toi and Tr>2 of the cor- 
responding Fermi functions. We recapitulate the full set 
of used fit parameters in Supplementary tables HI [Til Hill 
and [IV] 
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Gqpc (e 2 /h) 


Tfit 


Tdi = T D2 (K) 


VgSVg (fJ,V) 








0.043 


-36 


0.10 


0.08 


0.046 


-30 


0.20 


0.18 


0.050 


-31 


0.30 


0.28 


0.057 


-30 


0.40 


0.41 


0.062 


-30 


0.50 


0.49 


0.065 


-30 


0.60 


0.61 


0.063 


-31 


0.70 


0.71 


0.063 


-30 


0.81 


0.80 


0.057 


-30 


0.91 


0.89 


0.054 


-28 


1 


1 


0.049 


-36 



Table III: Summary of parameters used to fit the data 
shown in Supplementary figure l4b (cooldown 2). The 

conductance Gqpc is measured. The applied QPC voltage 
bias is here SVd = —36 //V. The lever arm is r\a = 0.059. 



SVoifiV) 


Tdi (K) 


T D2 (K) 


VaSVa (pV) 


-57 


0.089 


0.077 


-56 


-37 


0.073 


0.061 


-36 


-18 


0.055 


0.054 


-15 





0.044 


0.044 





9 


0.046 


0.043 


7 


19 


0.057 


0.052 


15 


28 


0.064 


0.055 


24 


37 


0.073 


0.060 


32 


48 


0.088 


0.070 


44 


58 


0.094 


0.080 


53 



Table IV: Summary of parameters used to fit the data 
shown in Supplementary figure [5p. (cooldown 3). Here 
SVd is the applied voltage bias and the parameter Tfu is set to 
the measured GQpch/e 2 — 0.5. The lever arm is r\c = 0.062. 



IV. SUPPLEMENTARY DISCUSSION 
REGARDING DEVIATIONS BETWEEN DATA 
AND FREE ID CHIRAL FERMIONS THEORY 

In this section, we discuss two possible explanations of 
the observed deviations between the data and the free 
ID chiral fermions predictions that can be seen in Fig- 
ures 3d, 4c and 4d. The first one relies on additional 
active quantum dot levels. The second one on energy 
redistribution with the inner edge channel. 

A. Effect of additional quantum dot levels 

In the presence of more than one active quantum dot 
level, the straightforward Equation 1 does not hold. In 
the sequential tunneling regime, the master equation ap- 
proach applies^. This can result in non-intuitive phe- 
nomena. In particular, it was shown that the gate volt- 
age position of a Coulomb peak can be shifted by about 
its width when the temperature changes^. We found a 
similar phenomenon could change the gate voltage sep- 
aration between the two dips observed in presence of a 
double step energy distribution. 
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Figure 6: Possible explanation with additional QD lev- 
els of observed small deviations on cooldown one. 

Data shown as symbols are those obtained in cooldown one 
and also displayed in Figure 4c. Dashed lines correspond to 
the slope predicted by rjc = 0.057 but shifted horizontally 
by ±7 /iV. Such shifts are compatible with the presence of a 
second active QD level located at e ~ — 10 fieV, below the 
first level. 

By solving the master equation with non-equilibrium 
energy distributions, using a slightly modified version of 
the code provided by^, and with a second level of en- 
ergy e nearby the first active level, we could reproduce 
the observed difference between rjcSVo and SVd shown 
in Figure 4d and Supplementary figure [6] Such second 
level for cooldown one is not visible in the top left in- 
set of Figure 2 but its existence is suggested by data 
taken in slightly different conditions (not shown) and by 
the observed T s > T D at T = 30 mK (see Figure 2). 
The presence of such a level can reduce |<5Vg| by a con- 
stant offset for le^Vol > |e| if asymmetrically coupled to 
the source and drain electrodes. Supplementary figure [S] 
shows tjgSVg extracted from cooldown one (data also 
shown in Figure 4c) together with dashed lines of slopes 
as expected from rjQ = 0.057 but offset by ±7 fiV which 
is compatible with a second level located at e w — 10 /J.eV, 
below the first level. 



B. Effect of energy redistribution with inner edge 
channel 

More recent data that we obtained using the same ex- 
perimental configuration show that an energy redistribu- 
tion exists between the outer and inner edge channels on 
length scales larger than 2 fim 17 . This effect, although 
probably small for the short 0.8 /im propagation length 
considered here, could explain the also small deviation 
observed between the free ID chiral fermions model and 
our data in Figures 3d and 4d. Another indication that 
energy relaxation along the edge is not fully negligible 
is the observed increase in the fit temperatures Tdi and 
Tdi with injected power (see Supplementary tables HI HT1 
lITTIandllVl). 
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